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A microscopic mean-field description of heavy ion collisions is performed in the framework of the 
time dependent Hartree-Fock theory using a Skyrme energy density functional. A good agreement 
with experiments is obtained on the position of the fusion barriers for various total masses and 
mass asymmetries. The excitation function of the ^^O+^'^^Pb is overestimated by about 16% above 
the barrier. Transfer below the barrier is studied in ^^0-|-^°*Pb central collisions. Heavier systems 



are considered to study their fusion hindrance. We also compute collision times of the U-|- 



U 



system. The latter has been used to produce super-strong electric fields and to test non perturbative 
quantum electrodynamics theory. Indeed, if the life time of such giant system is of the order of few 
10~^^ s, its electric field should induce spontaneous electron-positron pair emissions from vacuum. 
In our calculations, highest collision times are reached in the 238-[j^_238]|j j-gg^gi^jojj fQj- center of mass 
energies between 1000 and 1300 MeV. 



I. INTRODUCTION 

Description of nuclear reactions is very challenging, es- 
pecially at energies around the fusion barrier generated 
by the competition between Coulomb and nuclear inter- 
actions. It has been established that the structure of 
the collision partners may affect strongly the reaction 
mechanisms in this energy domain, as, for instance, the 
fusion cross-sections (for a review, see, e.g., Ref. 
In fact, a good understanding of these effects is cru- 
cial for the quest of superheavy elements (SHE) In- 
deed, different structures of the collision partners may 
change SHE production cross sections by several orders of 
magnitudes [3[. For instance, fusion hindrance has been 
observed for collisions of quasisymmetric heavy systems 



with typical charge products Zi x Z2 > 1600 [4j. Here, 
dynamical shape evolution of the dinuclear system and 
transfer between fragments are expected to play a sig- 
nificant role on this fusion hindrance. Investigations on 
collision dynamics between actinides are also of interest, 
in particular to study the physics of superstrong electric 
fields gig. 

It is therefore recommended to treat both structure of 
nuclei and their collision dynamics within the same for- 
malism. This is the case of fully microscopic approaches 
such as the time dependent Hartree-Fock (TDHF) theory 
proposed by Dirac in 1930 This time dependent ver- 
sion of the well known Hartree-Fock (HE) theory [l3| 
gives a self-consistent mean-field description of nuclear 
dynamics. The success of the first HE calculations based 
on Skyrme interactions [111, [l2| led to tremendous ac- 
tivities to describe nuclear structure within mean-field 
based approaches (see, e.g., Ref. [l^] for a review). Early 
TDHF calculations in nuclear physics [ij, [l^ have been 
done with the seek for a description of the dynamics of 
nuclei as good as their static properties [l^. At that 
time, theoreticians had to use various symmetries and 
simplified Skyrme interactions to reduce computational 



times. Now, recent increase of computational power al- 
lowed realistic TDHF calculations of nuclear collisions in 
3 dimensions with full Skyrme energy density functional 
(EDF) [13, [11. 

In this proceeding, we present a review of various 
TDHF studies of heavy ion collisions. In Section [Hi 
we recall the formalism and detail the calculation. In 
Section lllli we study fusion and transfer in light and 
medium heavy systems. In particular, a detailed study 
of the ^^0-|-^°^Pb reaction is presented. In Section [TVl 
we investigate the problem of fusion hindrance of heavy 
quasisymmetric systems. Then, we study the collision 
time of two ^asy nuclei in Section IVl before to conclude. 



II. THE TIME-DEPENDENT HARTREE-FOCK 
APPROACH 

In nuclear physics, the TDHF theory is applied with 
a Skyrme EDF modeling nuclear interactions between 
nucleons [11]. The EDF is the only phenomenological in- 
gredient of the model, as it has been adjusted on nuclear 
structure properties like infinite nuclear matter and radii 
and masses of few doubly magic nuclei [l^ . The main 
approximation of the theory is to constrain the many- 
body wave function to be an antisymmetrized indepen- 
dent particles state at any time. The latter ensures an 
exact treatment of the Pauli principle during time evo- 
lution. Though TDHF does not include two-body col- 
lision term, it is expected to treat correctly one-body 
dissipation which is known to drive low energy reaction 
mechanisms as Pauli blocking prevents nucleon-nucleon 
collisions. 

The main advantage of TDHF is that it treats static 
properties and dynamics of nuclei within the same for- 
malism and the same EDF. The initial state is obtained 
through static HE calculations which give a good approx- 
imation of nuclear binding energies and deformations. 
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Another important advantage of TDHF for near-barrier 
reaction studies is that it contains all types of couplings 
between the relative motion and internal degrees of free- 
dom whereas in coupled channels calculations one has to 
include them explicitly according to physical intuition. 
However, TDHF gives only classical trajectories for the 
time-evolution and expectation values of one-body ob- 
servables. In particular, it does not include tunneling of 
the many-body wave function. 

Inclusion of pairing correlations responsible for su- 
perfluidity in nuclei have been done recently within 
the time dependent Hartree-Fock-Bogolyubov theory to 
study pairing vibrations in nuclei [20| . However, realistic 
applications to heavy ion collisions are not yet achieved 
and are beyond the scope of this work. 



A. Formalism 

The TDHF equation can be written as a Liouville-Von 
Neumann equation 

^n^^p=[h[p],p] (1) 

where p is the one body density matrix associated to the 
total independent particles state with elements 

Ai+A-2 

p{Ysqys'q')= J2 ^.{^sq)v*{v's'q'). (2) 

1=1 

The sum runs over all occupied single particle wave func- 
tions Lpi and r, s and q denote the position, spin and 
isospin of the nucleon respectively. The Hartree-Fock sin- 
gle particle Hamiltonian h[p] is related to the EDF, noted 
E[p\, by its first derivative h[p\{Ysq, r's'q') = ^„(i.wlsol ■ 



B. Practical aspects 

A TDHF calculation of two colliding nuclei is per- 
formed assuming that the two collision partners are ini- 
tially at a distance Dq in their HF ground state. This 
distance has to be big enough to allow Coulomb excita- 
tion in the entrance channel (polarization, vibration, ro- 
tation...). We chose Dq = 44.8 fm in sections lllll and HVl 
and Do = 51.2 fm in sections fVl We assume that be- 
fore to reach this distance, the nuclei followed a Ruther- 
ford trajectory which determines their initial velocities. 
A Galilean transformation [13] is then applied on each 
nucleus in the first iteration. 

We use the tdhfSd code built by P. Bonche and 
coworkers with the SLyAd Skyrme parametrization p7j . 
This code has been extensi vely used to study heavy ion 
fusion [13, m m, m, m, . it has a plane of sym- 
metry (the collision plane). It uses the Skyrme energy 
functional expressed in Eq. (A. 2) of Ref. |28| where the 



tensor coupling between spin and gradient has been ne- 
glected. The lattice spacing is Ax = 0.8 fm and the time 
step is At 1.5 X 10"^"^ s. 



III. FUSION AND TRANSFER IN LIGHT AND 
MEDIUM-HEAVY SYSTEMS 

In this section, we focus on reactions with spherical 
nuclei. Studies of the effects of deformation on fusion 
within TDHF can be found in [H,!!!,!!!,!!^. 



A. Fusion barriers 

Fusion barriers are classically defined as the energy 
threshold above which fusion occurs for a head-on colli- 
sion. Experimentally, the average position of the barrier 
can be approximated by the centroid of the so-called bar- 
rier distribution [soj . To determine fusion barriers from 
TDHF, we consider head-on collisions at various ener- 
gies. The barrier of a system is then located between the 
highest energy for which there is no fusion and the lowest 
one for which fusion occurs. 

Figure [T] shows the density plot for a ^^0-|-^"*Pb 
central collision at a center of mass energy Ecm — 
74.44 MeV (top) and Ecm = 74.45 MeV (bottom). We 
see that the system separates into two fragments after 
a neck formation at the lowest energy. However, adding 
10 keV is enough to fuse. We deduce the fusion bar- 
rier V^°"^ = 74.445 ± 0.005 MeV for this system. This 
value is in excellen t ag reement with the experimental one 
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FIG. 1: Densities associated to a 0-f Pb central colli- 
sion at a center of mass energy E = 74.44 MeV (top) and 
74.45 MeV (bottom). The surface corresponds to an isoden- 
sity at half the saturation density. Density decreases from red 
to black. Each plot is separated by 135 fm/c. 

To see if these results are affected by dynamical effects, 
we compare them with a purely static approach. Fig- 
ure [5] shows the nucleus-nucleus potential of this system 
if we assume frozen HF densities at various internuclear 
distances R using the same EDF and the same numeri- 
cal approximations as in the TDHF calculations. Here, 
repulsion due to Pauli principle is neglected, which is a 
rather good approximation for distances down to the bar- 
rier as overlap of the nuclear densities are small. With 
this approach, the barrier is y^''"^^" — yg MeV [27j . 
This is shghtly higher than the TDHF value. This dif- 
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ference is due to dynamical effects sucli as transfer (see 
Section lill Bp that reduce tfie barrier. 
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FIG. 2: Nucleus- nucleus potential of "0-h^''*Pb assuming 
frozen HF densities of the nuclei at various internuclear dis- 
tances R. 



B. Transfer below the barrier in ^''0+^°*Pb 

It has been suggested that transfer may affect the bar- 
rier in the ^^0+^°®Pb system jsl]. Indeed, we see in 
Fig. [T] that the two fragments are linked by a neck and 
form a dinuclear system during ~ 400 fm/c. It is enough 
time for the nuclei to transfer nucleons, leading to a dy- 
namical evolution of the barrier. 

To get a deeper insight into the transfer preceding fu- 
sion in this reaction, we have plotted the evolution of the 
numbers of protons and neutrons of the small fragment as 
function of energy in Figure 21 Almost two protons and 
no neutron, in average, have been transfered from the 
^^O to the ^''^Pb. The two-proton transfer from the light 
to the heavy nucleus is then expected to be an important 
channel at the barrier. This is consistent with experi- 
mental observations of relatively high Carbon production 
cross sections, of the same order of the Nitrogen ones, in 
the exit channel of ^^O-h^^^Pb at the barrier [H [Hi]. 
Note that we have to distinguish between the transfer of 
typically one or two nucleons and more violent collisions 
like deep-inelastic reactions. In the latter, the width of 
the particle number distribution is known to be underes- 
timated with a single Slater determinant [36|. 



Figure [3] shows a comparison between experimental fu- 
sion barriers and those from TDHF calculations for sys- 
tems with various total masses and mass asymmetries. 
The lowest barrier is for ^°Ca-|-'"'Ca and the highest one 
for ^^Ti-|-^°*Pb. The agreement with experimental data 
is slightly better with TDHF than with the Bass bar- 
rier [34I. Remembering that these TDHF calculations 
have no adjustable parameter on reaction mechanisms, 
we conclude that one can use it for fusion barriers pre- 
dictions with confidence. 
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FIG. 3: Theoretical fusion barriers from TDHF calculations 
(circles) and the Bass barriers (stars) as function of the ex- 
perimental values (from barrier distributions centroid). 
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FIG. 4: Transfer below the barrier in ^^O-t-^^^^Pb. The aver- 
age number of protons and neutrons of the small fragment are 
plotted as function of energy. The arrow indicates the TDHF 
fusion barrier. 



C. Fusion cross section of ^'^0+2''*Pb 

We now focus on fusion cross sections. In TDHF, 
the latter can be approximated by <Jfus{EcA4) — T^b^ 
where be is the critical impact parameter associated to 
the initial Coulomb trajectrory below which fusion oc- 
curs [2^. Figure [5] shows the excitation function obtained 
for the ^^0-|-^°^Pb system in comparison to experimen- 
tal data [3l|. There is a good agreement above the bar- 
rier, though the fusion cross sections are overestimated 
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by about 16 %. The origin of this small discrepancy is 
nonetheless puzzling as the TDHF fusion barrier agrees 
perfectly with experiment. More precise investigations of 
fusion above the barrier with TDHF are then needed. 

Below the barrier, the TDHF fusion cross section van- 
ishes, following a classical behavior. This drawback of 
TDHF is well known and is due to the restriction to a 
single independent particles state. Indeed, to get a fusion 
probability between and 1, we need at least two Slater 
determinants: one describing the two well separated frag- 
ments after the collision when fusion does not occur and 
one describing the fused system. It is necessary to go 
beyond TDHF to treat a sum of Slater determinants and 
then to describe sub-barrier fusion due to quantum tun- 
neling. 



description of nuclear dynamics, should be appropriate 
to study fusion mechanisms in heavy (quasi)-symmetric 
systems. 
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FIG. 5: Experimental fusion excitation function (stars) for 
i6Q_^208p^j [Pi rpj^g jjjjgg denote the upper and 

lower limits for the fusion cross sections obtained from TDHF 
calculations. 



IV. FUSION HINDRANCE OF 
QUASISYMMETRIC HEAVY NUCLEI 

As we showed in the preceding section, TDHF system- 
atics are in good agreements with experimental fusion 
barriers for light or medium-heavy systems. In these 
cases, models where the nuclei are assumed to be frozen 
with their ground state density give already a good ap- 
proximation to fusion barriers. However, when collision 
partners get heavy and symmetric, typically when the 
product of their charges is such that Z\ Z2 > 1600, 
such frozen approaches significantly underestimate fu- 
sion barriers [4|. These systems need an extra energy 
to fuse compared to standard barrier height estimations. 
This phenomenon gave birth to the extra-push model [13] ■ 
One-body dissipation of the kinetic energy towards nucle- 
onic degrees of freedom is often invoked to be responsible 
for this fusion hindrance. Thus, TDHF theory, which in- 
corporates such one-body dissipation in the microscopic 
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FIG. 6: Fusion barriers obtained for various systems with the 
proximity model 38] (dotted line) , the extra-push model [s^] 
(dashed line) and TDHF (solid line). 



Let us now present a preliminary study of *®Ti-|-^°^Pb, 
56Fe+208pb, 9"Zr-|-i24Sn, ^^Zi+^^'^Xe and i24Sn-)-i36xe 
head-on collisions. Due to computational time limita- 
tions, we assume that the nuclei have fused if a com- 
pact system is formed for a minimum duration of 1200 
fm/c. Moreover, HF ground statse of ^^''Sn and ^^^Xe 
are slightly deformed. As a first approximation, we ne- 
glect possible effects of deformation on fusion and con- 
sider only one orientation to limit computational effort. 

Results are shown on Fig. [6l and compared with the 
extra-push model [STj and the proximity predictions of 
fusion barriers As we can see, TDHF results are 

in rather good agreement with the extra-push model. 
Both approaches predict significantly higher fusion bar- 
riers than the proximity model when mass symmetry in- 
creases. TDHF should then be useful tool to investigate 
dissipative mechanisms responsible for fusion hindrance. 

Fusion models usually rely on the reduction of the com- 
plex dynamical many-body problem into few relevant 
degrees of freedom, like neck parameters, mass asym- 
metry or distance between fragments. The TDHF the- 
ory, however, is free of such choice of collective vari- 
ables. Thus, the shape of the system can evolve freely in 
three dimensions. For example, we show on Fig. [7] snap- 
shots of the density profile in the reaction plane of the 
90Zr+i24Sn central coUision at Ecm = 235 MeV. This 
bombarding energy is 17 MeV higher than the barrier 
given by the Hartree-Fock frozen density approximation 

frozen ^ y^^y^ ^^j^^ TDHF fu- 

sion threshold (Vj^^^ = 238.75 MeV). The separation 
of the system in two fragments is thus entirely driven by 
the fusion hindrance mechanisms. 

We observe in Figure [7] the formation of a neck with a 
density close to saturation density. Also, this neck sur- 
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FIG. 7: Snapshots of the density profile in the reaction plane 
of the ^OZr+^^^Sn central collision at Ecm = 235 MeV. Hot 
colors are associated to densities close to the saturation den- 
sity po — 0.16 fm~'^. 

vives during approximately 1000 fm/c. This time is long 
enough to allow for a strong rearrangement of the sys- 
tem via nucleon transfer and dynamical shape evolution 
which modifies the fusion potential landscape. However, 
we also observe in Fig. [7] that this rearrangement is lim- 
ited as it does not allow for the two fragments to merge 
and form a single one, i.e., the system keeps clearly the 
shape of a dinuclear system. These observations are 
encouraging and studies of isospin equilibration, mass 
transfer and shape evolution during the contact configu- 
ration are ongoing. 



V. COLLISION TIME OF ACTINIDES 

Collisions of two Uranium or heavier nuclei are used, 
amongst other things, to investigate the physics of super- 
strong electric fields [(?]. Indeed, it has been proposed 
that a giant system with total charge Z > 173 may induce 
spontaneous electron-positron pair emissions from vac- 
uum by a fundamental quantum electrodynamics (QED) 
process Q. However, no experimental evidence of 
this process has been obtained so far [s^. The keypoint 
might be a good understanding of collision dynamics, in 
particular a reliabe prediction of collision times. In the 
238u_j_238y reaction, one might be able to see this pro- 
cess for sticking times of several 10 s as shown by re- 
cent theoretical calculations based on the time-dependent 
Dirac equation 7*|. 

Figure [8] shows that there is no pocket in the nucleus- 
nucleus potential of the 238u_|_238u gyg^gj-^^ agreement 
with Refs. [13, |4l|. However, dissipation mechanisms 
such as evolution of nuclear shapes may delay the sepa- 
ration of the system [i^l . Experimentally, collision times 
have been estimated from oscillations in (5— electron spec- 
tra [i^ and values of ~ 2 x 10^^^ s have been obtained in 
238u_|_238u above the barrier Q. Recently, delay times 
in this reaction was searched analyzing kinetic energy loss 
and mass transfer [i^. On the theoretical side, macro- 
scopic models have first been used (see, e.g., Ref. [i^l. 
However, the complexity of reaction mechanisms and the 



high number of degrees of freedom to be included in re- 
alistic calculations motivate the use of microscopic ap- 
proaches. First microscopic calculations of the collision 
of two 238u n^gigj have been performed recently thanks 
to the Quantum Molecular Dynamics (QMD) model [4]| . 
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FIG. 8: Symmetric fission valley of the 2asu_^2j8u ^^^^^ 
system. The Hartree-Fock (solid line) and Hartree-Fock- 
Bogoliubov (circles) adiabatic potential energies under a 
quadrupole constraint are plotted as function of the total 
quadrupole moment. 

Here, we report preliminary results on central collisions 
of two 238]j deformation axis parallel to each other 
and perpendicular to the collision axis. The collision time 
Tcoii between nuclei is defined, here, as the time dur- 
ing which the neck density exceeds jOo/10 = 0.016 fm^^. 
Figure M shows the evolution of Tcou as a function of 
Ecm- We see a rise and fall of Tcoii with a maximum of 
- 1000 fm/c at 1000 < Ecm < 1300 MeV, in a rather 
good agreement with the QMD calculations of Ref. [4l[ . 
Central collisions in this energy range are then expected 
to produce superstrong electric fields with lifetimes which 
fill the necessary condition to observe e^e~ pair excita- 
tions from QED vacuum. 



VI. CONCLUSIONS AND PERSPECTIVES 

We presented some applications of the TDHF theory 
to heavy ion collisions around the barrier. The only phe- 
nomenological input is the set of parameters of the en- 
ergy density functional which have not been adjusted on 
any reaction mechanism. Despite this, the agreement be- 
tween the TDHF calculations and the experimental fu- 
sion barriers is excellent for a wide range of projectiles 
and targets. However, these calculations overestimate 
the fusion cross sections for the system ^^0-|-^°^Pb above 
the barrier by about 16%. Theoretical investigations are 
needed to understand this (rather small) disagreement. 

Though the TDHF fusion cross section has a classical 
behavior, the quantum nature of the single-particle wave 
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FIG. 9: Collision time of two ^'^^U nuclei as function of center 
of mass energy. 



functions is well treated. Thus, we presented a study of 
transfer in the ^^O+^^'^Pb below the barrier. Two-proton 
transfer from the ^^O to the ^ospb nucleus is dominant at 
the barrier. This qualitative observation is in agreement 



with experiment. 

We are presently extending the fusion study to heavy 
quasi-symmetric systems to investigate the phenomena 
responsible for the fusion hindrance in such systems. Pre- 
liminary results are encouraging as they indicate that this 
fusion hindrance is well reproduced by TDHF. 

We finally presented calculations of collision times be- 
tween two 238^ j^^clei. Our conclusion is that the giant 
system might survive enough time (~ 1000 fm/c at an 
energy 1000 < Ecm < 1300 MeV) to test nonperturba- 
tive quantum electrodynamics with super-strong electric 
fields. 
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